TRIANGULAR EQUILIBRIUM POINTS IN THE GENERALIZED 
PHOTOGRAVITATIONAL RESTRICTED THREE BODY 
PROBLEM WITH POYNTING-ROBERTSON DRAG 

B.S.KUSHVAH 1 AND B.ISHWAR 2 
(l)JRF UNDER DST PROJECT (2) P. I. DST PROJECT 
B.R.A. BIHAR UNIVERSITY MUZAFFARPUR-842001, INDIA 
EMAIL: BSKUSH@HOTMAIL.COM 



Abstract. In this paper, we have found the equations of motion of Generalized 
Photogravitational restricted three body problem with Poynting-Robcrtson 
drag. The problem is generalized in the sense that smaller primary is sup- 
posed to be an oblate spheroid. The bigger primary is considered as radiating. 
The equations of motion are affected by radiation pressure force oblateness 
and P-R drag. We have located triangular equilibrium points in our problem. 
All classical results involving photogravitational and oblateness in restricted 
three body problem may be verified from this result. 
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1. Introduction 

In general three body problem, we study the motion of three finite bodies. The 
problem is restricted in the sense that one of the three masses is taken to be so 
small that the gravitational effect on the other masses by third mass is negligible. 
The small body is known as infinitesimal mass and remaining two massive bodies 
as finite masses or primaries. The classical restricted three body problem is gener- 
alized to include the force of radiation pressure, the Poynting- Robertson drag effect 
and oblateness effect. 

Poynting(1903) has stated that the particle such as small meteors or cosmic dust 
are comparably affected by gravitational and light radiation force, as they approach 
luminous celestial bodies. He also suggested that infinitesimal body in solar orbit 
suffers a gradual loss of angular momentum and ultimately spiral into the Sun. In a 
system of coordinates where the Sun is at rest, radiation scattered by infinitesimal 
mass in the direction of motion suffers a blue shift and in the opposite direction 
it is red shifted. This gives rise to net drag force which opposes the direction 
of motion. The proper relativistic treatment of this problem was formulated by 

Robertson(I937) who showed that to first order in — the radiation pressure force 
is given by 
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(1) 




where F v = 



P 16ttR' 2 psc 



3Lm 



denotes the measure of the radiation pressure force, R the 



position vector of P with respect to radiation sours S, V the corresponding velocity 
vector and c the velocity of light. In the expression of F p , L is luminosity of the 
radiating body, while m, p and s are the mass, density and cross section of the 
particle respectively. 

The first term in equation (0 express the radiation pressure. The second term 
represents the Doppler shift of the incident radiation and the third term is due to 
the absorption and subsequent re-emission of part of the incident radiation. These 
last two terms taken together are the Poynting- Robertson effect. The Poyinting- 
Robertson effect will operate to sweep small particles of the solar system into the 
Sun at cosmically rapid rate. 

Stanley P., Wyatt J. R. and Fred L. Whipple (1950) have shown that the P-R 
effect has been of very little significance galacticaly.Its importance is restricted to 
small particles orbiting in the vicinity of individual stars. The importance of the 
radiation influence on celestial bodies has been recognized by many scientists such 
as Kozai (1961), McCraken and Alexander (1968), Ferrz-Mello (1972), Simmons et 
o/.(1985), Vkrouhlixky (1993, 1994), Murray CD. (1994), Ragos O.and Zafiropoulos 
F.A. (1995). 

Colombo et al. (1996) analyzed the stability of the equilibrium points in the 
presence of radiation pressure which include the Poynting-Robertson drag term. 
They showed that the points were unstable to such a drag force. Chernikove (1970) 
has dealt with the Sun-Planet-Particle model. He concludes that the P-R effect 
renders unstable those liberation points known to be conditionally stable in the 
classical case. Schuerman D.W. (1980) has studied the triangular points. He has 
shown that these points are unstable on time scale long compared the period of 
revolution of the two massive bodies. 

This problem has an interesting application for artificial satellite and future 
space colonization. It has been suggested that the classical triangular points of the 
Sun-Jupiter or Sun-Earth system would be convenient sites to locate future space 
colonies. 

In this paper we consider Sun-Planet-Particle model with Sun as a radiating 
body, planet as an oblate spheroid. We have found that the coordinates of equi- 
librium points are the functions of mass reduction factor qi, P-R drag Wi and 
coefficients of oblateness Ai . All classical results in restricted three body problem 
involving radiation and oblateness may be deduced from this result. 

2. Equations of Motion 

We suppose mi, the mass of more massive radiating primary and mi the mass 
of smaller primary which is an oblate spheroid. Let these primaries revolve around 
their centre of mass in circular orbits in the plane of motion which coincides to the 
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equatorial plane of mi. 

We consider the barycentric rotating co-ordinate systemOXYZ relative to iner- 
tial system with angular velocity u> and common Z -axis. We have taken line joining 
the primaries as X -axis. OX and OY in the equatorial plane of mi and OZ coin- 
ciding with the polar axis of mi. Let r e andr p be the equatorial and polar radii of 
mi and r be the distance between primaries. Let infinitesimal mass m be placed at 
a point P{X, Z, Y").Thcn potential at P due to m\ and mi is 

k 2 mm\q\ k 2 mmi k 2 mmiAi 

( ' ~ n ~ 2 2r| 

r 2 — r 2 

where k 2 is the Gaussian constant of gravitational and Ai = e 5r2 p ,oblatcncss coef- 
ficient of mi, r\ and ri are the distances from m\ to P and mi toP respectively. 
We take units such that sum of the masses and distance between primaries as unity. 
The unit of time i.e. time period of mi about mi will consists of 2tt units such 
that k = 1. Then perturbed mean motion of the primaries is n 2 = 1 + ^j 2 -. Let 
u = — ^ — then — ^ — = 1 — a with mi > m 2 where u is mass parameter. Then 
(HI and (El and use the same technique as Chernikove (1970) and Schuerman 
(1980). In dimensionless coordinate system, the dimensionless velocity of light will 
be given by Cd = c, which depends on the physical masses of the primaries and the 
distance between them. 

For inertial reference system the total acceleration acting on P is given by Coriolis 
relation 

(1 - /i)rl [Lfi 3 \iAifi 



2uj x v + lu x (lu x r) 



4 2 r| 



(1 - /i)(l - q\) f rl (r'l + uj x rl).rlrl r\ + lu x rl 



ri c d r 2 c d 

where 

r"i = + yj, v = xi + yj, a — xi + yj , lu — nk, f{ — (x + fi)i + yj, 

r 2 = (x + \i - 1)1 + yj, r\ = (x + p,) 2 + y 2 , r\ = (x + // - l) 2 + y 2 . 

Substituting all these values in above relation and comparing the components of 
i and j ,we get the equations of motion of the infinitesimal mass particle in xy-pl&ne. 



U x =x-2ny = r 2 x ^ - + ^ ^ + ^ - ^ 3 M 2 (a; + /i-l) 



/on Wl f (* T + M) r/ s • -i 

(3) T < 2 — [{x + (J,)x + yy\ + x - ny 



2 



T j •• , o • 2 (1 - fJ.)giy ny 3 nA 2 y 
U v = y + Inx = n y 5 5 — - — g — 

Wi ( y 

(4) r <^ — [(x + fx)x + yy] + y + n(x + fi) 

rf ( rj 

where W\ = ilziflilzlll 77,2 = 1 + ^A 2 ,q = 1 — & is a mass reduction factor 
expressed in terms of the particle radius a, density p radiation pressure efficiency 
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factor x ( m C.G.S. system) q = 1 — 5,6 **° — \- The assumption g=constant is 
equivalent to neglecting fluctuations in the beam of solar radiation and the effect 
of the planets shadow. Obviously q < 1. Equations of motion (OJ) and (2} can be 
written as 

a; - 2m/ = + F x , y + 2nx = + F y 



where 



U 1 = ^(x 2 + y 2 ) + i l^ + » + ^ 
2 r\ r 2 2r 2 



/r\ J — i Wi f (x + u) r/ , . .. 

(5) F x = T < j — [{x + fi)x + yy\ + x - ny 



Wi ( y 

(6) and F v = ^- < \{x + fj,)x + yy] + y + n(x + fi) 

r( {r( 

F x ,F y are the partial derivatives of drag function with respective to x and y respec- 
tively, which are purely functions of the particle's position and velocity. 
Now multiplying equations ([5J by 2i;,(|fjJ by 2y and adding, we get, 

2xx + 2yy^2(x-^-+y-^-j + 2(xF x +yF y ) 

which can be written as ^ = — 2{xF x + yF y ), where C = 2Ui — x 2 — y 2 . The 
quantity C is Jacobi Integral. The zero velocity curves are given by C = 2U\(x,y) 

3. Location of Triangular Equilibrium Points 

For the triangular equilibrium points y ^ 0, U x = U y = then from equations 
(0 and (i) 

(7) nx + -r«y = o 

'1 '2 z '2 '1 



( 8 ) " y la — 3 — ra(x + /x) = 



(l-fi)qxy \iy 3 /iA 2 y Wi 
2 r| ri 

Multiplying equations (EJ by y, ( IB} by (a; + /i) then subtracting, we get 



{ n2 -4 (l + ^ ) } W = nWl - 



In the case of photogravitational restricted three body problem we have rx — 
q 1 / 3 = S (say) and r 2o = q 2 ^ 3 — 1-So we suppose that due to P-R drag and 
oblateness, perturbation in rx and r 2o are ex and e 2 respectively where e/s are 
very small, 

(9) rx = gj /3 (l + ex) and r 2 = 1 + e 2 

Putting these values in equation ( [7J and ( [SJ and neglecting higher order terms 
of small quantities, we get 

nWx 

e2V ~ 3M1 + |A 2 ) 
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nWi(l-fA 2 ) nW 1 A, 
we may write as e 2 y = e 2 y , i-e. £2 = 5 > £ i = ~7T; ^ 7T Hence, 

we get 

f nWi A 2 \ . x * 5 

ri =5U- -r; r - f,r 2 = 1 + (1 - -A 2 ) 

I 6(1 - /i)y 2 J 3^y 2 



Since x + /1 = — — ^ — and y 2 = r\ — (x + /i) 
(10) x = x \ 1 



nWM(l- M )(l-|A 2 )+/x(l-^)ir] ^2 
3/i(l - ^)2/oa;o 2 a; 

_ f n^ 2 [2 M -l-Ml-M^)f +7(1- M )^] ^ 2 (l-f)A 2 ^ 1/2 

(ii)y - yo r smi-^o 3 ^ 

(%o,yo) ar e coordinates of £4,^5 in the photogravitational restricted three body 
problem where 



52 54 V /2 & 

Xo = — -H, Vo = ±0 1 1 - -r I , = 9i 



1/3 



Equations ( ITHjl and ( [TTJ) are valid for Wy < 1, j4 2 <C 1. 

4. Conclusion 

The equilibrium points L4, and L5 are given by equations ( 110(1 and ( 111(1 . Case 
1. We see that when P-R effect not included i.e.,Wi = then we get x — xo — -^A 2 
and y = yo{l — ^"(1 — ir)y^} which are the L4 and £5 points in the case of 
generalized photogravitational restricted three body problem. Case 2. When 
A 2 = then n = 1, i.e., small primary is spherically symmetric. We have 



3^(l-^)x yo /' °\ 6/i(l-/i)y§ 



2 



this result coincides with Schuerman(1980). Case 3. When A 2 = 0, Wi = 0,then 
we have for q\ = 1 = S, x = | — (1, y = These are the coordinates of classical 
restricted three body problem. 

Finally we conclude that position of triangular equilibrium points are affected 
by mass reduction factor, P-R drag and oblateness coefficient. 

Acknowledgment : 

We are thankful to D.S.T. Government of India New Delhi for sanctioning a 
project DST/MS/140/2K on this topic. 

References 

[1] Beauge C. and Ferraz-Mello S.(1994): Capture in exterior mean-motion resonances due to 

Poynting-Robertson drag. Icarus 110, 239-260. 
[2] Chernikov Yu.A.(1970): The photogravitational restricted three body problem. Sov. Astr. 

AJ. 14( 1), 176-181. 

[3] Colombo G., D.A. Lautman and Shapiru 1.1.(1996): The Earth's dust belt:Fact or fiction? 

2. Gravitational focusing and Jacobi capture. J. Geophys. Res. 71, 5705-5717. 
[4] Ishwar B.(1997): Nonlinear stability in the generalized restricted three body problem. Celst. 

Mech. and Dyn. Astron. 65, 253-289. 



6 B.S. KUSHVAH AND B. ISHWAR 

[5] Kalacka J. (1992): Poynring-Robertson effect 1. Equations of motion. Earth Moon and Planets, 
Kluwer Academic Publishers, Printed in the Neitherland 59, 41-59. 

[6] Liou J.C., Zook and Jackson(1995): Radiation pressure, Poynting-Robertson drag and solar 
wind drag in the restricted three body problem. Icarus 116, 186-201. 

[7] McCuskey S.W.(1963):Itroduction to celestial Mechanics. Addison Wesley Publishing com- 
pany Inc. New York. 

[8] Murray C.D.(1994):Dynamical effect of drag in the circular restricted three body problem 1. 

Location and stability of the Lagrangian equilibrium points. Icarus 112, 465-484. 
[9] Poynting J. H. (1903): Radiation in the solar system: Its effect on temperature and its pressure 

on small bodies. Phil. Trans. Roy. Soc. London, 202-523. 
[10] Ragos O. and Zafiropoulos F. A. (1995): A numerical study of the influence of the Poynting- 
Robertson effect on the equilibrium points of the photogravitational restricted three body 

problem l.Coplanar case. Astro. Astrophys. 300, 568-578. 
[11] Robertson H.P.(1903): Dynamical effects of radiation in the solar system.M.on. Not. Roy. 

Astron. Soc. 97, 423-438. 
[12] Sahoo S.K.and Ishwar B.(2000): Stability of collinear equilibrium points in the generalised 

photogravitational elliptical restricted three body problem. Bull. Astro. Soc. India 28, 576- 

586. 

[13] Shuerman D.W.(1980): The restricted three body problem including radiation pressure. As- 
trophys. J. 238(1), 337-342. 

[14] Simmon J.F.L., Mc Donald A.J. Brown J.C.(1985):The restricted three body problem with 
radiation pressure. Celst. Mech. 35, 145-187. 

[15] Singh J. and Ishwar B.(1999):Stability of triangular points in generalised photogravitational 
restricted three body problem. Bull. Astron. Soc. India 27, 415-424. 

[16] Wyatt S.P.and Whipple F.L.(1950): The Poynting-Robertson effect on meteor orbits. Amer- 
ican Astron. Soc. APJ. Ill, 134-141. 



